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Abstract: 

The complete set of electroweak 0{a) corrections to the Drell-Yan-like production of 
W bosons is calculated and compared to an approximation provided by the leading term 
of an expansion about the W-resonance pole. All relevant formulae are listed explicitly, 
and particular attention is paid to issues of gauge invariance and the instability of the 
W bosons. A detailed discussion of numerical results underlines the phenomenological 
importance of the electroweak corrections to W-boson production at the Tevatron and at 
the LHC While the pole expansion yields a good description of resonance observables, 
it is not sufficient for the high-energy tail of transverse-momentum distributions, relevant 
for new-physics searches. 
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1 Introduction 



The Drell-Yan-like production of W bosons represents one of the cleanest processes with 
a large cross section at the Tevatron and at the LHC This reaction is not only well suited 
for a precise determination of the W-boson mass M-w, it also yields valuable information 
on the parton structure of the proton. Specifically, the accuracy of ~ 15-20 MeV [Q] in the 
Mw measurement envisaged at the LHC will improve upon the precision of ~ 30 MeV to be 
achieved at LEP2 and Tevatron Run II [Q, and thus competes with the precision of the 
Mw measurement expected at a future e"^e~ collider Concerning quark distributions, 
precise measurements of rapidity distributions provide information over a wide range in 
a; 0; a measurement of the d/u ratio would, in particular, be complementary to HERA 
results. The more direct determination of parton-parton luminosities instead of single 
parton distributions is even more precise extracting the corresponding luminosities 
from Drell-Yan-like processes allows one to predict related qq processes at the per-cent 
level. 

Owing to the high experimental precision outlined above, the predictions for the pro- 
cesses pp/pp W Ivi should match per-cent accuracy; for specific observables the 
required theoretical accuracy is even higher. To this end, radiative corrections have to 
be included. In particular, it is important to treat final-state radiation carefully, since 
photon emission from the final-state lepton significantly changes the lepton momentum, 
which is used in the determination of the W-boson mass. A first step to include elec- 
troweak corrections was already made in Ref. |^, where effects of final-state radiation in 
the W-boson decay stage were taken into account. Those effects lead to a shift in the value 
of Mw of the order of 50-150 MeV. The approximation of Ref. [0 was improved much 
later in Ref. where the electroweak 0{a) corrections to resonant W-boson produc- 
tion P were discussed for W production at the Tevatron in detail. The 0{a) corrections 
are of the order of the known next-to-next-to-leading order (NNLO) QCD corrections ||1( 



to Drell-Yan-like processes. In 0{a'^) only a study of two-photon radiation exists |jTT 
while the virtual counterpart is completely unknown. Discussions of QCD corrections to 
Drell-Yan-like processes can be found in Refs. [1^,|13| and references therein. 



In this paper we present the complete calculation of the electroweak 0{a) corrections, 
including non-resonant contributions. In particular, we compare the full 0{a) correction 
to a pole approximation (similar to the one used in Refs. that is based on the correc- 

tion to the production of resonant W bosons. All relevant formulae are listed explicitly. 
Moreover, a discussion of numerical results is presented for the Tevatron (Run II) and for 
the LHC. Partial results of this analysis have already been presented in the LHC workshop 
report [0 . 



The paper is organized as follows. In Sect. ^ we set our conventions and provide 
analytical results for the parton-level subprocess. In particular, we describe the calculation 
of the complete and the "pole-approximated" 0{a) corrections. Different methods for 
treating the infrared and coUinear singularities are presented and compared with each 
other. The hadronic cross section is discussed in Sect. ^ In Section ^ we present numerical 
results for W-boson production at the Tevatron and at the LHC. Our conclusions are given 
in Sect. ||. Finally, the appendices provide some supplementary formulae. 
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Figure 1: Lowest-order diagram for ud — > W+ — > uil^. 

2 The parton process ud — > W+ — > ^z^^(+7) 

2.1 Conventions and lowest-order cross section 

We consider the parton process 

uiPu,'T'u) + diPd,rd) ui{kn,-) +l^iki,Ti) [+7(A;, A)], (2.1) 

where u and d generically denote the hght up- and down-type quarks, m = u, c and d = d,s. 
The lepton / represents / = e, /i, r. The momenta and hehcities of the corresponding 
particles are given in brackets. The Mandelstam variables are defined by 

s = {Pu+ Pdf, i = {pd - kif, u = {pu- hf, Suii = {kn + kif. (2.2) 

We neglect the fermion masses m^j, m^, m/ whenever possible, i.e. we keep these masses 
only as regulators in the logarithmic mass singularities originating from coUinear photon 
emission or exchange. Obviously, we have s = Si,^i for the non-radiative process ud vil^ ■ 
In lowest order only the Feynman diagram shown in Fig. [T| contributes to the scattering 
amplitude, and the Born amplitude reads 



with an obvious notation for the Dirac spinors Vd-, etc., and the left-handed chirality 
projector = |(1 — 75). The electric unit charge is denoted by e, the weak mixing angle 
is fixed by the ratio = 1 — = M-^/M| of the W- and Z-boson masses Mw and Mz, 
and Vud is the CKM matrix element for the ud transition. 

Strictly speaking, Eq. ( |2.3| ) already goes beyond lowest order, since the W-boson width 
Fw(s) results from the Dyson summation of all insertions of the (imaginary part of the) 
W self-energy. Defining the mass Mw and the width Fw of the W boson in the on-shell 
scheme (see e.g. Ref. [|15|), the Dyson summation directly leads to a running width, i.e. 

rw(S)U = Fw^. (2.4) 

On the other hand, a description of the resonance by an expansion about the complex 
pole in the complex s plane corresponds to a constant width, i.e. 

rw(S)|,„„,t = Tw- (2.5) 
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In lowest order these two parametrizations of the resonance region are fully equivalent, 
but the corresponding values of the line-shape parameters Mw and Fw differ in higher 
orders i,|T6|,0, 



const 



w 



T.2 , 

J- W + ■ ■ 
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(2.6) 



Since Mwlmn— Afwiconst ~ 26 MeV, it is necessary to state explicitly which parametrization 
is used in a precision determination of the W-boson mass from the W line-shape. 

The differential lowest-order cross section is easily obtained by squaring the lowest- 
order matrix element A^o of 
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|S-M2, + iMwrw(s)|2 



(2.7) 



where the explicit factor 1/12 results from the average over the quark spins and colours, 
and fl is the solid angle of the outgoing in the parton centre-of-mass (CM) frame. The 
electromagnetic coupling a = e^/(47r) can be set to different values according to different 
input-parameter schemes. It can be directly identified with the fine-structure constant 
a(0) or the running electromagnetic coupling a(Q^) at a high-energy scale Q. For instance, 
it is possible to make use of the value of a(M|) that is obtained by analyzing |TB| the 
experimental ratio R = cr(e"'"e~ — *• hadrons)/(e"'"e~ — * fi^fi~). These choices are called 
a{0)-sche'me and a{M^) -scheme, respectively, in the following. Another value for a can 
be deduced from the Fermi constant G^, yielding ac^ = V^G^M^s'^/tt; this choice is 
referred to as G ^j^- scheme. The differences between these schemes will become apparent in 
the discussion of the corresponding 0{a) corrections. 



2.2 Virtual corrections 



The virtual one-loop corrections comprise contributions of the transverse part of the 
W self-energy SJ^, corrections to the two Wdu and Wvil vertices, box diagrams, and 
counterterms. The explicit expression for SJ^ (in the 't Hooft-Feynman gauge) can, e.g., 
be found in Ref. 
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The diagrams for the vertex and box corrections, which are shown in 
Fig. ^, were calculated using standard methods. The Feynman diagrams and amplitudes 
were generated with FeynArts [|19|. The subsequent algebraic reduction of the one- 



loop tensor integrals to scalar integrals was performed with FeynCalc and the scalar 



integrals were evaluated using the methods and results of Ref. The algebraic part was 
checked numerically by a completely independent calculation, in which the amplitudes are 
expressed in terms of tensor coefficients using Mathematica, and the tensor reduction is 
done numerically. UV divergences are treated in dimensional regularization, and the IR 
singularity is regularized by an infinitesimal photon mass m^. The actual calculation is 
performed in 't Hooft-Feynman gauge using the on-shell renormalization scheme described 
in Ref. [|l^, where, in particular, all renormalization constants used in this paper can be 
found. As an additional check we have repeated the calculation within the background- 



field formalism ||23| and found perfect agreement. In the following we sketch the structure 
of the virtual corrections and emphasize those points that are relevant for the treatment 
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Figure 2: Diagrams for vertex and box corrections. 



of the resonance and for the change from one input-parameter scheme to another. The 
complete expressions for the vertex and box corrections are provided in App. 

While the self-energy and vertex corrections are proportional to the Dirac structure 
appearing in the lowest-order matrix element Aio, the calculation of the box diagrams 
leads to additional combinations of Dirac chains. However, since the box diagrams are 
UV-finite, the four-dimensionality of space-time can be used to reduce all Dirac structures 
to the one of Aio, see App. |A.2|. In summary, the complete one-loop amplitude Aii can 



be expressed in terms of a correction factor 5™^ times the lowest-order matrix element. 

Ml = 6^'''Mo. (2.8) 

Thus, in 0{a) the squared matrix element reads 

\Mo + Mi\'' = (l + 2Re{5^''^*})|A^or + ..., (2-9) 

so that the Breit-Wigner factors are completely contained in the lowest-order factor | A^oP- 
Note that the Dyson-summed imaginary part of the W self-energy, which appears as rw(s) 
in A4o, is not double-counted, since only the real part of enters Re{5"''*} in 
The correction factor 5"''* is decomposed into four different parts, 

5"''* = Swwis) + Swduis) + SwuA-^) + ^hcis, t), (2.10) 



^Note also that only the imaginary parts induced by the fermion loops, which lead to the physical 
decay width in the resonance propagator, are resummed, while all other parts of are treated in 0{a). 
In this way the differences between the on-shell renormalization of the W-boson mass (see e.g. Ref. |l5| ) 
and other variants based on the complex pole position in the inverse W propagator are of 0{a^), and thus 
beyond the accuracy of our calculation. In particular, the problems with the resummation of loops 
containing photon exchange are avoided. 



according to the splitting into self-energy, vertex, and box diagrams. 
The W self-energy correction reads 

= _ ^T(g)-jM^ _ ^z^^ (2.11) 

where the explicit expression for the unrenormalized self-energy is given in Eq. (B.4) 
of Ref. |jl5|. In the on-shell renormalization scheme the renormalization constants for the 
W-boson mass and field, ^Zw, are directly related to SJ^. 

The vertex corrections are given by 

6wdu{s) = Fwdu{s) + 6^duy Swuii{s) = Fwuilis) + S^^^i, (2.12) 

where the explicit expression for the form factor Fwff'{s) is given in App. |A.1| . The 
counterterm for the Wff vertex depends on the input-parameter scheme. In the 

a(0)-scheme (i.e. the usual on-shell scheme), it is given byQ 



+ h6Zw + 6Zf'^ + 6Z^''^). (2.13) 



The wave-function renormalization constants SZ'^'^ and 6Z-^''^ are obtained from the (left- 
handed part of) the fermion self-energies, and the renormalization of the weak mixing 
angle, i.e. 5sw, is connected to the mass renormalization of the gauge-boson masses. The 
charge renormalization constant 6Ze contains logarithms of the light-fermion masses, in- 
ducing large corrections proportional to aln{m'j/s), which are related to the running of 
the electromagnetic coupling a((5^) from Q = to a high-energy scale. In order to ren- 
der these quark- mass logarithms meaningful, it is necessary to adjust these masses to the 
asymptotic tail of the hadronic contribution to the vacuum polarization of the photon. 
Using a(M|), as defined in Ref. |jl8|], as input this adjustment is implicitly incorporated, 
and the counterterm reads 



^H>//'L(M|) = ff Lio) - ^^"(^z), (2-14) 



where 

-)2\ -nAA /r,\ o„rTT^A/^2 



Aa(Q^) = (0) - Re{n;- (Q^)}, (2.15) 

with IIj^^ denoting the photonic vacuum polarization induced by all fermions other 
than the top quark (see also Ref. [|15|)- In contrast to the a(0)-scheme the counter- 



term 5vi///'|Q,(jvf2) does not involve light quark masses, since all corrections of the form 
a"ln"(mj/s) are absorbed in the lowest-order cross section parametrized by q;(M|) = 
q;(0)/[1 — Aa(M|)]. In the G^-scheme, the transition from a(0) to is ruled by the 
quantity Ar |l^,^, which is deduced from muon decay, 

= V2G,M^sl ^ ^^^^^^ ^ ^ ^^^3^_ ^2.16) 

TC 



^We consistently set the CKM matrix to the unit matrix in the correction factor 5'^"^, since mixing 
effects in the 0{a) corrections are neghgible. This means that the CKM matrix appears only in the global 
factor iKidP to the C'(a)-corrected parton cross section. 
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Therefore, the counterterm S^^j-, reads 

Since Aa(M|) is exphcitly contained in Ar, the large fermion-mass logarithms are also re- 
summed in the G^-scheme. Moreover, the lowest-order cross section in G^-parametrization 
absorbs large universal corrections induced by the p-parameter. 

The box correction (5box(S, t) is the only virtual correction that depends also on the 
scattering angle, i.e. on the variables t and u = —s — t. The explicit expression for 
Shox{s,t) is given in App. |A.2 . 



Despite the separation of the resonance pole (s — M-^)~^ from the correction factor 
^virt (|2.8|) , 6^^^^ still contains logarithms ln(s — + that are singular on resonance. 
Since these singularities would be cured by a Dyson summation of the W self-energy inside 
the loop diagrams, we substitute 

ln(s-M^ + ie) ln(s - + iMwFw) (2.18) 

with a fixed width everywhere, independent of the use of a fixed or running width in 
lowest order. In principle, also a running width could be used on the r.h.s., but the 
difference to the fixed width is of two-loop order, and thus beyond the accuracy of our 
calculation. The substitution ( |2.18| ) does not disturb the gauge-invariance properties of 



the one-loop amplitude Aii, i.e. its gauge-parameter independence and the validity of 
SU(2)xU(l) Ward identities. The reason is that the sum of all terms in A4i proportional 
to ln(s — M-^ + ie) separately fulfills the (algebraic) relations arising from gauge invariance, 
since this logarithm is the only resonant (and thus a unique) term in the relative correction 

^virt_ 



2.3 Virtual correction in pole approximation 

If one is only interested in the production of (nearly) resonant W bosons, the electroweak 
corrections can be approximated by an expansion about the resonance pole, which 
is located in the complex s plane at — iMy^T^^ up to higher-order terms. The ap- 
proximation of taking into account only the leading term of this expansion is called pole 
approximation (PA) and should not be confused with the on-shell approximation for the 
W bosons. In contrast to the PA, the on-shell approximation, where the W bosons are 
assumed to be stable, does not provide a description of the W line-shape. In the following 
we construct a PA for the virtual correction factor in the same way as a double-pole 



approximation was constructed in Ref. for the more complicated case of W-pair pro- 
duction in e"'"e~ annihilation, e"'"e~ WW 4fermions. In this formulation the PA is 
only applied to the virtual corrections, while the real corrections are based on the full 
photon-emission matrix element, as described in the next section. In principle, it is also 
possible to construct a PA for the real corrections, as for instance described in Ref. |^ 



for ee —>■ WW.0 However, we prefer to make use of the full matrix elements, as it was also 
done in the variant of the PA used in Ref. |Q] for W production in hadronic collisions. 

•^More details about the pole expansions that are used in practice are reviewed in Ref. [E9|. 



6 




on-shell production on-shell decay 

Figure 3: Generic diagram for factorizable corrections. 



In the PA the virtual corrections to wrf W — > uil'^ can be classified into two cat- 
egories. The first category comprises the corrections to the production and the decay of 
an on-shell W boson. Owing to the independence of these subprocesses, these corrections 
are called factorizable. All contributing Feynman graphs are of the generic form shown 
in Fig. 1^. By definition, the factorizable corrections receive only contributions from the 
W self-energy and the W f f vertex corrections. The corresponding correction factor 5l[^^^ 
is obtained from 5ww{.s) and 5vk//'(s) by setting s = and Fw = 0. Since we have 
Sww{M^)\ry^=o = in the complete on-shell renormalization scheme, we get 



rvirt 
"fact 



>Wdu 



rw=0' 



:2.i9) 



Since the vertex corrections for on-shell W bosons correspond to physical S-matrix ele- 
ments, both contributions to the factorizable corrections are gauge- invariant. Note that 
6flll is a constant factor, neither depending on the scattering energy nor on the scattering 
angle. Moreover, S'^^l contains IR singularities originating from the logarithms of ( |2.18D ; 
these terms are connected to photon emission from on-shell W bosons and are regularized 
by the infinitesimal photon mass m^. 



The second category of corrections in the PA are called non-factorizable pU|,pT| and 
comprise all remaining resonant contributions, i.e. all terms in 5^"^* that are non-vanishing 
for s — > Mw and Fw 0, 



^nonfact('^) ^) 



rvirt 
''fact- 



(2.20) 



Here (5"''*|s^jv/2^,rw^o results from the full off-shell correction 5"''* upon taking the asymp- 
totic limits (s — Myj) and Fw — * while keeping the ratio (,§ — M^)/(MwFw) fixed. 
As can be shown by simple power counting [^|, only loop diagrams with an internal pho- 
ton contribute. The relevant diagrams for ud ^ W ^ vil"^ are shown in Fig. ^ Since 
also box diagrams are involved, production and decay do not proceed independently, and 
the terminology non-factorizable is justified. The limit s — > for the pole expansion 
has to be defined carefully, because s is not the only kinematical variable. The variables t 
and u, which are related to the scattering angle 6 of the parton CM frame, range within 
—s<t,u<0 and are related to s by s + 1 + -u = 0. Therefore, changing s while keeping t 
and u fixed is inconsistent in general. We circumvent this problem by taking s 
fixed scattering angle 6, resulting in the replacements 



M2,for 



t tr 



t 



u 



Ur 



U ■ 



-M^cos^f. 



(2.21) 
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Figure 4: Diagrams for non-factorizable corrections (y? is the would-be Goldstone partner 
to the W boson, and u±, u^ denote Faddeev-Popov ghosts). 



The actual calculation of ^nonfact is performed as described in Ref. in detail. The final 
result is 

Oact(S, i) = -f (-2 + Li, (l + ^) - Qu Li2 (l + ^ 



+ 2 In 

where 



- iMwF 



w 



m^Mw 

dt 



g^ln --^ -Q„ln 



^2.22) 



Li2(x) = - / — In(l-t), |arc(l -a;)| < vr, (2.23) 
JO t 

is the usual dilogarithm, and Qf denotes the electric charge of the fermion /. In ( |2.22| ) 
we made use of Qi = Qd — Qu- We note that ^nonfact is, by definition, a gauge-invariant 
quantity. It does not involve mass singularities of the fermions, but it is IR-singular. 
More precisely, the IR-singular term proportional to aln(m^) exactly compensates the 
artificially created IR singularity in 6^111 by setting s = and Fw — > there, and these 
a In(m^) terms of ^^act replaced by the correct logarithms a ln(s — + iMwFw) after 
C'nfact is added. 

In summary, the PA for the virtual correction factor reads 

< = '^faJt* + C„\act(^;i) (2.24) 



with Sj^ll and (5nonfact given in ( |2.19| ) and ( ^.221 ), respectively. The uncertainty induced by 



omitting non-resonant corrections can be estimated naively by 

< - ~ f ^ ^, ln( ^ ) , In^ ( ^ ) ^ (2.25) 
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Figure 5: Diagrams for real-photon emission. 



where the factor Fw/Mw resuhs from the neglect of non- resonant contributions near 
resonance and the logarithms account for typical enhancements away from resonance. 

Finally, we note that we have analytically compared the PA for the virtual corrections 
worked out in this section with the results presented in Ref. 0. Apart from non-resonant 
contributions, which go beyond the validity of the PA, both results agree. 



2.4 Real-photon emission 

Real-photonic corrections are induced by the diagrams shown in Fig. ^. In the calcula- 
tion of the corresponding amplitudes it is mandatory to respect the Ward identity for 
the external photon, i.e. electromagnetic current conservation. Writing the amplitudes as 
Ai^{X) = e*^{k, X)T^ with e* denoting the polarization vector of the outgoing photon, this 
Ward identity reads k^T'^ = 0. If the W width is zero, this identity is trivially fulfilled. 
This remains even true for a constant width, since the W-boson mass appears only in the 
W propagator denominators, i.e. the substitution — ^ — iM^Fw is a consistent 
reparametrization of the amplitude in this case. However, if a running W width is intro- 
duced naively, i.e. in the W propagators only, the Ward identity is violated. The Ward 
identity can be restored by taking into account the part of the fermion-loop correction 
to the jWW vertex that corresponds to the fermion loops in the W self-energy leading 
to the width in the propagator In Ref. |Q it was shown that this modification 



simply amounts to the multiplication of the 'jWW vertex by the factor 

WU = l + i^ (2-26) 



if the photon is on shell {k = 0). By construction [|r7|,R2|,p3], the width Fw appearing in 



f-ywwlrun ^^^^ ^ propagator is the lowest-order width, since it results from 

the imaginary part of the W self-energy at one loop. Since, however, for the special case 
of an on-shell photon the relevant imaginary parts are completely parametrized by the 
ratio Fw/iWw, the Ward identity is fulfilled for any numerical value of Fw- Therefore, 
we are allowed to use a value for Fw that includes also QCD and electroweak radiative 
corrections. For later convenience, we define 

flWwlconst = 1- (2.27) 

In addition to the U(l) Ward identity for the external photon, an SU(2) Ward iden- 
tity ||17| , ^ becomes relevant for effectively longitudinally polarized W bosons in W7 



"^There is a misprint in Eq. (D.45) in Ref. the factor 2 in front of the Qi'Q f'[{f;i) — > 
Qi'Q f[{i,t) — > {i',u)] and QiQ f'[{.f,t) {f',u)] terms needs to be removed. 
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production at high energies, where the lepton current [u^^'jfj^uj^vi] becomes proportional to 
the W momentum few = kn + h. As can be checked easily, this Ward identity is maintained 
in the above treatment of the W width as well. 

The helicity amplitudes tVI^"^'*^' (A) for the radiative process ud — >• can be written 
in a very compact way using the Weyl-van der Waerden spinor formalism. Adopting the 
conventions of Ref . 133] , we obtain 



M: 



+ 



QuiknkiY 



[sy^i - M^j + lM■sNT^^{sy^l)\{puk)(J)dk) 



s-M^ + iMwrw(s) 



Qi{PuPd)* ^ (Qd - Qu){puki) 

{kkn){kki) {pdk){kkn) 



+ 



f-yWwiQu - Qd){knkl)*{Puk)* 



M: 



(-1) 



+ 



+ 



{{PuhYf 
1 



Qdiknk) 



[s,,i - + iMwTwM]{puk)*{pdky 

QuiiPuPd) , {Qd- Qu){Pdkn) 



s-M^ + iMwrw(s) 



{kkn)*{kkiY {pukY{kkiY 



f-yWwiQd - Qu){knkl){pdk) 



(2 28) 

[s-M2, + iMwrw(s)][s,,;-M2, + iMwrw(s,,0](P«^)*i' 

The amplitudes for the other helicity channels vanish for massless fermions. The spinor 
products are defined by 



{pq) 



e^^PAQB 



'^VPo% 



9p 
2 



'^^ COS -f sm 



2 



'''' cos — sin — 



2 



(2.29) 



where pa, Qa are the associated momentum spinors for the light-like momenta 

P^ = Po(l) sin 9p cos (pp, sin 9p sin 0p, cos 9p) , 
q'^ = Po(l) sin6'q cos0g, sin6'g sin0g, cos6'g). 



(2.30) 



The contribution of the radiative process to the parton cross section is given by 



= — — / dr^ y\M: 

^ 12 2s J 



XA)P 



(2.31) 



where the phase-space integral is defined by 

, d^K f d^ki r d^k 
dr. ' ' 



i27TY2kn,oJ (27r)32A;i,o V (27r)32A;o 



{2ny6{pu+Pd-kn-k-k). (2.32) 



2.5 Treatment of soft and collinear singularities 



The phase-space integral ( |2.31|) diverges in the soft (fco 0) and collinear {puk,pdk, kik — > 
0) regions logarithmically if the photon and fermion masses are set to zero, as done in 
( p.28| ). For the treatment of the soft and collinear singularities we applied three different 
methods, the results of which are in good numerical agreement. In the following we briefly 
sketch these approaches. 
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2.5.1 IR phase-space slicing and effective collinear factors 

Firstly, we made use of the variant of phase-space slicing that is described in Ref. ||35| , 
where the soft-photon region is excluded in the integral ( |2.31| ) but the regions of photon 
emission collinear to the fermions are included. 

In the soft-photon region < < AE <^ \/I the bremsstrahlung cross section 
factorizes into the lowest-order cross section and a universal eikonal factor that depends 



on the photon momentum k (see, e.g., Ref. ||T5[). Integration over k in the partonic CM 
frame yields a simple correction factor (5soft to the partonic Born cross section dao. 



■'soft 




2AE\ 



+ 2 In 



+ 



\ J \ s 
\ J \ s 



2 In 



2 In In 



+ 21n 

+ 2 In 
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/2AE\ 
\ ) 
/2AE\ 
\ J 
f2AE\ 



+ 



IT 



2 In In 



-u 



Li2 

Li2(- 



2 



+ - W 



+ 



IT 



U 

i 
i 
u 



(2.33) 



The factor 5soit can be added directly to the virtual correction factor 2 Re{(5"''*} defined in 
(p. 91). It can be checked easily that the photon mass cancels in the sum 2 Re{5"''*}+5soft- 
The remaining phase-space integration in (|2.31|) with k^ > AE still contains the 
collinear singularities in the regions in which {puk), (pdk), or (kik) is small. In these 
regions, however, the asymptotic behaviour of the differential cross section (including its 
dependence on the fermion masses) has a well-known form. The singular terms are uni- 
versal and factorize from dcro- A simple approach to include the collinear regions consists 
in a suitable modification of |7Vl^p, which was calculated for vanishing fermion masses. 
More precisely, is multiplied by an effective collinear factor that is equal to 1 up to 

terms of 0{m'j/s) (/ = /, u, d) outside the collinear regions, but replaces the poles in (puk), 
{pdk), and (kik) by the correctly mass- regularized behaviour. Explicitly, the substitution 
reads 

Y: |-M;-^-^KA)r ^ E ft\^u,Xu,E^,e^,)fif\md,Xd,Ed,ed,) 

A=±l K„,Kd,K;=±l 

X fj;''\mi,xi,Ei,ei^) |A^^"-"''^''--'^'-'(A)p. (2.34) 

A=±l 

The functions and f^^^/^'^^ describe collinear photon emission with and without 

spin flip of the radiating fermion, respectively, 

ft^'-\mf,Xf,Ef,ej 



AE} sin\efj2) 



AEjsm^{ef^/2) + mjJ ' 



f^-^'^\mj,Xf,Ef,9f,) 



X 



f 



Am)E) sin2(^^^/2) 



X, 



T'^Xf + 2[AE'jsm.\efj2) + 



m 



212' 



Xf = ^. (2.35) 
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where Ef is the fermion energy and 6f^ = Z(kj,k) is the angle of the photon emission 
from f = u,d, I. 



2.5.2 IR and collinear phase-space slicing 



Instead of using effective collinear factors, alternatively we have also applied phase-space 
slicing to the collinear singularities, i.e. the collinear regions are now excluded by the 
angular cuts Of^ < A9 <^ 1 in the integral ( ^.31[ ). The IR region is treated as previously, 
leading to the same correction factor 4oft as given in ( |2.33| ). 

In the collinear cones the photon emission angles 6f^ can be integrated out by making 
use of the factorization property of the squared photon-emission matrix elements with the 
radiator functions j^™'/^"^)^ described in the previous section. The resulting contribution 
to the bremsstrahlung cross section has the form of a convolution of the lowest-order cross 
section. 



C"coll 



C"coll,n + CTcolLd + CTcoUJ, 



27T Jo 



2n Jo 



1-2A_B/Vl 




dz 



+ {l-z)cro{ki) 



(2.36) 

q = u,d, (2.37) 



with the splitting function 



l + z' 



(2.38) 



(2.39) 



For initial-state radiation the respective quark momentum pq is reduced by the factor z so 
that the partonic CM frame for the hard scattering receives a boost, while this is not the 
case for final-state radiation. Note that for final-state radiation, i.e. in cxcoii,/, the lepton 
momentum in the final state is zki, although ki is relevant for the lowest-order cross section 
in the hard scattering. Of course, if photons collinear to the lepton / are not separated 
the z integration can be carried out explicitly. It can be checked easily that in this case 
all logarithms of the lepton mass m/ cancel in the sum of virtual and real corrections. 



2.5.3 Subtraction method 



Finally, we applied the subtraction method presented in Ref. [^, where the so-called 
"dipole formalism", originally introduced by Catani and Seymour [^] within massless 
QCD, was applied to photon radiation and generalized to massive fermions. The general 
idea of a subtraction method is to subtract and to add a simple auxiliary function from 
the singular integrand. This auxiliary function has to be chosen such that it cancels all 
singularities of the original integrand so that the phase-space integration of the difference 
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can be performed numerically. Moreover, the auxiliary function has to be simple enough 
so that it can be integrated over the singular regions analytically, when the subtracted 
contribution is added again. 

The dipole subtraction function consists of contributions labelled by all ordered pairs 
of charged external particles, one of which is called emitter, the other one spectator. For 
ud — > z/^/^ we, thus, have six different emitter/spectator cases //': ud, du, ul, lu, dl, Id. 
The subtraction function that is subtracted from J2x l-^-yWl'^ is given by 



\Msnh\^ = QuQde^ gud{Pu,Pd,k)\MBorn{XudPu,Pd,kl,ud)\^ 

+ gdu{Pd,Pu, k)\MBoTn{Pu, XudPd, kl,du)\^ 

- QuQie^[gui{Pu, k, k) + giu{kuPu, k)\ \MBoTn{xuiPu,Pd, h^ui)? 
+ QdQie^[gdi{Pd, ki, k) + gid{kuPd, k)\ \MBorniPu,XdiPd, h^di)?, 

with the functions 

gud{Pu,Pd, k) 



(2.40) 



1 



gdu{pd,Pu,k) = 

gql{Pq,kl,k) = 

giuiki,Pq,k) = 

and the auxiliary variables 
PuPd - Puk - Pdk 



Xud 
1 

{Pdk)xud 
1 

{Pqk)Xql 
1 

{kik)xqi 



1 Xud 

2 

1 Xud 

2 



2 Xql Zql 

2 

2 Xql Zql 



Xud 
Xud 
-I- Xql 
-1- Zql 



Xud 



X, 



PuPd 



ql 



Pqkl +Pqk - kik 
Pqkl + Pqk 



Zql 



q = u,d, 



Pqkl 
Pqkl + Pqk ' 



[2.41] 



[2.42) 



For the evaluation of |A^subP in ( |2.4CI| ) the lepton momenta kijfi still have to be specified. 
They are given by 



k^^^, = A{pu,Pd)^k^, = k{pd,Pu)\k\, 

with the Lorentz transformation matrix 



(P + PY{P + P)y 2P^P, 



+ 



Km = K + /c^ - (1 - XqM. (2.43) 



p'' = ki" + k';:, p^ = XudPt + f2- 



P2 + PP ■ P2 ' 

(2.44) 

The modified lepton momenta kijfi still obey the on-shell condition kf ffi = 0, and the 
same is true for the corresponding neutrino momenta that result from momentum con- 
servation. It is straightforward to check that all coUinear and soft singularities cancel in 
X)a I-^7('^)P~ l-^subP so that this difference can be integrated numerically over the entire 
phase space (|2.32| ). If phase-space cuts are applied to the lepton momentum, these cuts 
directly affect ki in \M..^{X)\'^, but in |A1subP they have to be applied to hjfi. The 
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singularities nevertheless properly cancel in this case, since the momenta kijf are defined 
in such a way that they asymptotically approach ki in the singular regions. 

The contribution of |A1subP, which has been subtracted by hand, has to be added 
again. This is done after the singular degrees of freedom in the phase space ( p.32| ) are 
integrated out analytically, keeping an infinitesimal photon mass and small fermion 
masses mj as regulators |36|. The resulting contribution is split into two parts: one that 
factorizes from the lowest-order cross section o"o and another part that has the form of a 
convolution integral over (Tq with reduced CM energy. The first part is given by 



a 



0"sub,l — T;"^ QuQd 

Ztt 



+ QiQ 

with the auxiliary function 



C{-i, ml ml) 



2^ 
1 7r2 



QiQu 



1 vr^ 



C{-u,mi,mJ - - - — 



(2.45) 



C{r, mi,m2) = 2 In 



mim2\ ^^(m^ 



-) 



In - 

V r 



+ 21n 

\ r 



In^ 



-In 



fmim2\ 

\^)- 
(2.46) 

The IR and fermion-mass singularities contained in d(Tsub,i exactly cancel the ones of the 
virtual corrections. The second integrated subtraction contribution is given byQ 



crsub,2(p«,Pd) = QuQd— / dx <^ [Qud{s,x)] ao{xpu,Pd) 

Ztt Jo 

-QiQu^ [ du [ dx [giu{-u,x)]^{xpu,Pd 
ZTT J-s Jo au 



S = (pu+Pd)^ 

U={pd-kn)'^ 



+ QiQd 



a 



27r J- 



dt / dx 



d(To, 



Gid{-t,x) —^{Pu,xpd) 
-'+ dt 



t — {pu ^n)'^ 



+ ^ [ dx (1 - x) lQl(ro{xpu,Pd) + Qd(^o{Pu, xpd) 

Ztt Jo I tu^-tu 



where the usual 



prescription. 



dx 



fix] 



g{x) 



dx/(x) [^(x) -^(1)] 



(2.47) 



(2.48) 



is applied to the integration kernels 

Ogq'{,r,x) = Pff{x) 
Glg{r,x) = Pff{x) 



In 



In 



^x(l — x)m'l 



(2.49) 



^Note that we did not literally follow the formulae of Ref. 
spin- flip parts. 



2(1 -x) 
but rearranged some terms in the 
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In ( p. 47 ) we indicated explicitly how the Mandelstam variable r has to be chosen in terms 
of the momenta in the evaluation of the part containing [Qffi{r,x)]^. Note, however, that 
in ( p.47| ) the variable s that is implicitly used in the calculation of (Jo(. . .) is reduced to 
2xpuPd- 

In summary, within the subtraction approach the real correction reads 



0-^ 



/ dr, 

12 2s 



EI-^7(A)P 

. A 



\M 



sub 



sub,l 



sub,2- 



(2.50) 



It should be realized that in Usub,! and o"sub,2 the full photonic phase space is integrated 
over. This does, however, not restrict the subtraction approach to observables that are fully 
inclusive with respect to emitted photons, but rather to observables that are inclusive with 
respect to photons that are soft or coUinear to any charged external fermion (see discussions 
in Sect. 6.2 of Ref. and Sect. 7 of Ref. ||3^). In this context, the implementation of 



phase-space cuts, which define the observable, in Eq. (|2.50|) is crucial. The phase-space 
cuts have to be applied in the subtraction parts Csub,!, crsub,25 and jA^subP, which are all 
defined on the 2^2 kinematics, exactly in the same way, in order to guarantee that 
these parts consistently cancel each other. On the other hand, the observable-defining 
cuts directly enter the 2 — > 3 part J2\ l-^-yWl"^- Consistency requires that these cuts treat 
soft or collinear photons inclusively, i.e. that photons coUinear to charged fermions have 
to be associated with these fermions; otherwise |A^subP cannot locally cancel the collinear 
singularities in J2\ l-^-yWl"^- 



3 The hadronic processes pp, pp — ^ W+ — > ^^^z(+7) 

The proton-(anti-)proton cross section a is obtained from the parton cross sections a^'^i'^^) 
by convolution with the corresponding parton distribution functions ^'1,2(2;); 

(^(■5) = / dxi / dx2qi{xi)q2{x2)cr^'^'''^\xiPg^,X2Pg2). (3.1) 
^ "'0 Jo 



Q1Q2 



In the sum Y.qiq2 the quark pairs qiq2 run over all possible combinations ud and du of 
up-type quarks m = u, c and down-type quarks ci = d, s. The squared CM energy s of the 
pp (pp) system is related to the squared parton CM energy s by s = X1X2S. 

The (9(a)-corrected parton cross section a^'i^'''^') contains mass singularities of the form 
a\n{mq), which are due to collinear photon radiation off the initial-state quarks. In com- 
plete analogy to the MS factorization scheme for next-to-leading order QCD corrections, 
we absorb these collinear singularities into the quark distributions. This is achieved by 
replacing q{x) in ( p.l|) according to 

g(x)-g(x,M2)- r-g(-,M2)f Q4p^^(;^)|ln(^)-21n(l-z)-l|l , (3.2) 

where M is the factorization scale. This replacement defines the same finite parts in the 
0{a) correction as the usual MS factorization in L)-dimensional regularization for exactly 
massless partons, where the In(mq) terms appear as 1/(1^ — 4) poles. In ( p.2| ) we have 
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preferred to exclude the soft-photon pole by using the [...]+ prescription. This procedure 
is fully equivalent to the application of a soft-photon cutoff, as used in Ref. [Q. 

The absorption of the collinear singularities of 0{a) into quark distributions, as a 
matter of fact, requires also the inclusion of the corresponding 0{a) corrections into the 
DGLAP evolution of these distributions and into their fit to experimental data. At present, 
this full incorporation of 0{a) effects in the determination of the quark distributions has 
not yet been performed. However, an approximate inclusion of the 0{a) corrections to 
the DGLAP evolution shows that the impact of these corrections on the quark distri- 
butions is well below 1%, at least in the x range that is relevant for W-boson production 
at the Tevatron and the LHC Therefore, the neglect of these corrections to the parton 
distributions is justified for the following numerical study. 



4 Numerical results 
4.1 Input parameters 

For the numerical evaluation we used the following set of parameters. 



a 


= 1/137.0359895, 


a(M|) 


= 1/128.887, 


= 1.16637 X 10-5 GeV"^ 


Mw 


= 80.35 GeV, 


Mz 


= 91.1867 GeV, 


Mh = 150 GeV, 


Fw 


= 2.08699 . . . GeV, 




= 0.119, 




rUe 


= 0.51099907 MeV, 




= 105.658389 MeV, 


rur = 1.77705 GeV, 


ruu 


= 4.85 MeV, 


rric 


= 1.55 GeV, 


mt = 174.17 GeV, 


nid 


= 4.85 MeV, 




= 150 MeV, 


mb = 4.5 GeV, 


iKdl 


= 0.975, 


\Vns\ 


= 0.222, 




iKdl 


= 0.222, 


\Vcs\ 


= 0.975, 


(4.i; 


which 


is consistent with experimental data 39| . Except for the CKM matrix elements 



this input is identical with the one used in the LEP2 Monte Carlo workshop report ||2^ 
on precision calculations for LEP2. 

We recall that the above set of data is overcomplete, but all the numbers are needed for 
the evaluation of the 0{a) corrected cross section in the different input schemes described 



in Sect. |2.1| . The masses of the light quarks are only relevant for the evaluation of the 
charge renormalization constant 6Ze, which drops out in the a{Mz)- and G^-schemes; 
the values for these masses are adjusted to reproduce the hadronic contribution to the 
photonic vacuum polarization The value for the W-boson decay width Fw is the 

0{a)- and 0(as)-corrected SM prediction in the scheme. Note that ag only enters the 
calculation of Fw in the results presented here, since we do not consider QCD corrections 
to the scattering processes. If not stated otherwise, the presented cross sections and 
distributions are calculated for a fixed W-boson width. We consistently take the CTEQ4L 
quark distributions for the evaluation of the pp and pp cross sections. If not stated 



otherwise the factorization scale M is set to the W-boson mass Mw, which is the standard 
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Figure 6: Total parton cross section a in parametrization and relative corrections 6 for 
different parametrizations; the respective PAs for 5 are shown for comparison. 



v^/GeV 


40 


80 


120 


200 


500 


1000 


2000 


ao/pb 


2.646 


7991.4 


8.906 


1.388 


0.165 


0.0396 


0.00979 


5/% 


0.7 


2.42 


-12.9 


-3.3 


12 


19 


23 


Spa/% 


0.0 


2.40 


-12.3 


-0.7 


18 


31 


43 



Table 1: Total lowest-order parton cross section ctq in parametrization and correspond- 
ing relative correction 6, exact and in PA. 



choice for resonant W production. A thorough investigation of the full QCD and QED 
factorization scale dependence would require the inclusion of the relevant QCD corrections 
and QED-improved parton densities, which is beyond the scope of this paper. Instead we 
at least study the dependence of the C(a)-corrected hadron cross sections on the QED 
factorization scale M using parton densities without QED evolution. 



4.2 Results for the parton process 

Figure |^ shows the total partonic cross section a and the corresponding relative correction 
S for intermediate energies. Note that the total cross section (including its correction) 
is the same for all final-state leptons / = e, /z, r in the limit of vanishing lepton masses. 
As expected, the parametrization of the Born cross section minimizes the correction 
at low energies, since the universal corrections induced by the running of a and by the 
p parameter are absorbed in the lowest-order cross section. Moreover, the naive error 
estimate ( |2.25| ) for the PA turns out to be realistic. The PA describes the correction in 
the resonance region within a few 0.1%. Table |I| contains some results on the partonic 
cross section and its correction up to energies in the TeV range. Far above resonance the 
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ffrun(Mw) I 1 1 1 1 1 1 r 

ffconst(Mw,const) _ ^ J ^ (^_^^) 

1.02 

1 

0.98 

Mw — Mw, const 

0-96 - Mw = Mw,™ - 

I I I I I I I 

60 65 70 75 80 85 90 95 100 

Vs/GeV 

Figure 7: Ratio of ud i/^/i"''(+7) cross sections in the G^-scheme evaluated with running 
and constant widths in the W propagator. 




PA cannot describe the exact correction anymore, since non-resonant corrections become 
more and more important. The difference between the full 0{a) correction and the PA 
is mainly due to (negative) Sudakov logarithms of the form aln^(s/M^), which are not 
contained in the PA, but in the full 0{a) correction. The large positive correction at high 
energies is induced by the radiative return to the W resonance via photon emission from 
the initial-state quarks and from the W boson, i.e. the reaction proceeds as W7 production 
with subsequent W decay in this case. 

Finally, we investigate the relation between the parametrizations of the W resonance 
by a constant or running width. As explained in Sect. at tree level the transition from 
a running to a constant width in the cross section is equivalent to the change ( |2.6| ) in 
Mw and Here we check whether, or to which accuracy, this statement remains valid 
in the presence of the 0{a) corrections. Figure |^ shows the ratio (3"run/o'const of the cross 
sections for ud — > z/^yU"^(+7), evaluated with a running width (o"run) and with a constant 
width ((Tconst)- Switching from one parametrization to the other without correcting Mw 
and Pw changes the cross section at the level of a few per cent near the resonance, as 
expected. However, adjusting additionally Mw and Pw according to ( |2.6| ), absorbs this 
change in the cross section up to a difference that is below 0.1%. This means that also 
in the presence of the 0{a) corrections the results of a W-mass or width determination 
can be transformed from one parametrization to the other via ( |2.6D , without repeating the 
analysis in the other parametrization. 

The results of this section have been obtained with the three variants described in 
Sect. ^]5| for combining virtual and real corrections. For slicing parameters 2AE/\^ and 
A6 in the range 10^^-10^^ the results agree within the integration errors, which are 
significantly below the accuracy of the numbers in Table |1|. 
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4.3 Results for pp W+ at the LHC 

We first consider W production at the LHC, i.e. we assume a pp initial state with a CM 
energy of ^/s = 14 TeV. For the experimental identification of the process pp W+ 
h'il~^{+j) we take the set of phase space cuts 

PT,/>25GeV, |^T>25GeV, h\ < 1-2, (4.2) 

where pt,z and rji are the transverse momentum and the rapidity of the charged lepton 
respectively, and = Pt,ui is the missing transverse momentum carried away by 
the neutrino. Note that these cuts are not "coUinear-safe" with respect to the lepton 
momentum, so that observables in general receive corrections that involve large lepton- 
mass logarithms of the form Q;ln(m;/Mw). This is due to the fact that photons within 
a small collinear cone around the charged lepton momentum are not treated inclusively, 
i.e. the cuts assume a perfect isolation of photons from the charged lepton. While this is 
(more or less) achievable for muon final states, it is not realistic for electrons. In order to 
be closer to the experimental situation for electrons, we additionally consider the following 
photon recombination procedure: 

1. Photons with a rapidity \ri^\ > 2.5, which are close to the beams, are treated as 
invisible, i.e. they are considered as part of the proton remnant. 



2. If the photon survived the first step, and if the resolution Ri^ = ^ {rji — rj^Y + (f)}^ 
is smaller than 0.1 (with 0;^ denoting the angle between lepton and photon in the 
transverse plane), then the photon is recombined with the charged lepton, i.e. the 
momentum of the photon and / are added and associated with the momentum of /, 
and the photon is discarded. 

3. Finally, all events are discarded in which the resulting momentum of the charged 
lepton does not pass the cuts given in 



While the electroweak corrections differ for final-state electrons and muons without photon 
recombination, called "bare" leptons in the following, the corrections become universal in 
the presence of photon recombination, since the lepton-mass logarithms cancel in this 
case, in accordance with the KLN theorem |^T|. From the explanations at the end of 



Sect. p.5.3| it is clear that the described subtraction approach is not applicable without 



photon recombination, i.e. for bare leptons, while the slicing variants can be used with 
and without photon recombination. 

Figures ^ and ^ show the distributions in the transverse-momentum px,/ and in the 
transverse- invariant-mass M^^^^i in pp — * W"^ vil^ {+'-)) for the LHC energy, together 
with the corresponding relative electroweak corrections 5. The transverse-invariant-mass is 



defined by M^^vii = y2pT,/|^T(l — cos0;y;i), where (pi/^i is the angle between the lepton and 
the missing momentum in the transverse plane. The distributions show the well-known 
kinks at px.z ~ Mw/2 and M^ j^^i ^ Mw, which are used in the W-mass determination. 
Near these kinks the correction 6 reaches the order of 10-20% for bare leptons, where the 
larger corrections occur in the electron case, because the logarithm aln(m;/Mw) is larger 
in this case. Since these enhanced corrections originate from collinear final-state radiation. 
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Figure 8: Transverse-momentum distribution (da/dpT^i) and relative corrections 6. 
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Figure 9: Transverse-invariant-mass distribution {da /dM^^i^^i) and relative corrections 5. 
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pp z///+(+7) at = 14TeV 



PT,//GeV 


25-00 


50-cx) 


100-cx) 


200-cx) 


500-00 


lOOO-oo 


cTo/pb 


1933.5(3) 


11.50(1) 


0.8198(4) 


0.1015(1) 


0.005277(1) 


0.0003019(1) 




-5.5(1) 


-9.2(1) 


-11.8(1) 


-16.5(1) 


-27.5(2) 


-39.1(1) 


<^e+i/e,PA/% 


-5.4(1) 


-8.0(1) 


-7.3(1) 


-7.8(1) 


-10.3(2) 


-13.0(1) 




-2.9(1) 


-4.9(1) 


-8.5(1) 


-13.1(1) 


-23.4(1) 


-34.5(1) 


'^^t+i/^,PA/% 


-2.8(1) 


-3.5(1) 


-4.0(1) 


-4.4(1) 


-6.2(1) 


-8.5(1) 


Srec/% 


-1.8(1) 


-2.7(1) 


-6.2(1) 


-10.2(1) 


-19.6(1) 


-29.6(1) 




-1.8(1) 


-1.5(1) 


-1.6(1) 


-1.6(1) 


-2.4(1) 


-3.6(1) 


A(2±iMw)/% 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.2 


A(Mt,.,z)/% 


-0.0 


0.0 


0.1 


0.2 


0.5 
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Table 2: Integrated lowest-order pp cross sections ctq for different ranges in px.i and corre- 
sponding relative corrections 6, exact and in PA. The QED scale uncertainty is illustrated 
by A(M), as described in the text. 



they are negative for higher pt,« and redistribute events to lower transverse momenta. The 
correction 6 is reduced to a few per cent after photon recombination, which eliminates the 
artificial lepton-mass logarithms. 

Moreover, Figures ^ and ^ demonstrate the reliability of the PA for transverse lepton 
momenta pt,i ^ Mw/2, where resonant W bosons dominate. The PA curves are included 
in the plots as thin double-dashed lines close to the corresponding full corrections. Only 
for transverse momenta above the resonance region a systematic difference between the 
PA and the full result starts to become apparent. 

Transverse lepton momenta pt,z above the resonance kink need to be considered for the 
determination of the W decay width Fw and, in the high-energy tail, for searches for new 
physics, such as new W gauge bosons. Table § shows the high-pT,/ contributions to the 
total cross section defined by different ranges in p^^i. Although the cross section rapidly 
decreases for high pt,^ sizeable event numbers can be expected for px,/ values in the range 
100-1000 GeV at the LHC. The results of Table | reveal that the PA is not apphcable for 
very large p^j, where the W boson is far off shell. This fact was also to be expected from 
the results at the parton level discussed in the previous section. As mentioned previously, 
the breakdown of the PA is mainly due to the missing Sudakov logarithms, which are 
independent of the lepton species. This explains why the PA breaks down both for electrons 
and muons in the same way, with or without photon recombination. Moreover, these large 
negative corrections of up to ~ 30% for p^ i values between 500-1000 GeV do not become 
small after photon recombination. Note also that the large positive correction to the 
total parton cross section at high scattering energies (see Table |l]), which is mainly due 
to resonant W'j production in the forward and backward directions, is widely suppressed 
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at the hadron level. The reason for this suppression is the cut on the transverse lepton 
momentum which effectively demands large W momenta. 

Table |^ additionally shows the uncertainty of the (!?(a)-corrected cross sections that 
is due to the dependence on the QED factorization scale M, as described in Sect. ^. The 
quantity 

A{M) = 5{M)-5 (4.3) 

quantifies the difference of the relative correction S{M), evaluated for the scale M, to 
its value 6 for the default choice M = Mw. This difference A(M) applies to each row 
in Table since the M-dependent part is the same for each version of 6,„ shown in the 
table. We have chosen the values M = 2^^Mw and M = M^^^^i for illustration, the 
former representing a simple rescaling of our default choice, the latter a scale choice that 
might be more appropriate for large transverse momenta. Although the parton densities 
have not been corrected for the QED factorization (see Sect. |), the QED scale dependence 
turns out to be negligible compared to the other theoretical uncertainties and the expected 
experimental accuracy. Using QED-improved parton densities the QED factorization scale 
dependence should be reduced further. 

The results without photon recombination have been obtained with the two slicing 
variants of Sect. |2.5| with the same values for the cut parameters as at the parton level. The 
results with photon recombination have also been checked with the subtraction approach. 
For the distributions shown in Figs. § and ^, the results from the different methods agree 
at the level of the line width in the plots.0 The results of Table § have been obtained 
with the various methods with an agreement within 1-2 standard deviations, apart from a 
few exceptions due to fluctuations (signalled by a large x^/d-o.f.). Apart from those, 
all numbers without photon recombination result from the "effective coUinear factor" 
approach, while the ones with photon recombination are obtained with subtraction. 

4.4 Results for pp — > W""" — > vil^ at the Tevatron 

Now we consider W production at the Tevatron, i.e. at a pp collider with a CM energy 
of a/s = 2 TeV. We again use the phase space of (f4.2|) and the photon recombination 
procedure^ of the previous section. 

The corrections to the transverse-momentum and transverse-invariant-mass distribu- 
tions at intermediate transverse momenta, displayed in Figs. ^ and |rT|, are of the same 
size as the ones discussed for the LHC in the previous section. They also show the same 
qualitative features under the changes in the lepton mass and after photon recombination. 
Moreover, the results for the PA are again not distinguishable from the full correction in 

^Of course, the statistical error of the histogram bins is somewhat larger than the line width, at least for 
large transverse momenta and invariant masses. This means that the differences between the subtraction 
and slicing variants are even smaller than the statistical error. However, this simply results from the 
fact that in all cases the same set of random numbers was used, which introduces correlations among the 
statistical errors in the contributions that are not affected by subtraction or slicing. 

^More realistic electron and muon identification procedures for the Tevatron have been used in Ref. § , 
where electroweak corrections to the transverse momentum and invariant-mass distributions have been 
evaluated in another variant of the pole approximation. 
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Figure 10: Transverse-momentum distribution {da/dpT,i) and relative corrections 6. 
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Figure 11: Transverse-invariant-mass distribution {da / dM^^^y^i) and relative corrections 5. 
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pp z///+(+7) at = 2 TeV 



PT,//GeV 


25-00 


50-cx) 


75-00 


100-cx) 


200-00 


300-00 


cTo/pb 


407.03(5) 


2.481(1) 


0.3991(1) 


0.1305(1) 


0.006020(2) 


0.0004821(1) 




-5.4(1) 


-9.7(1) 


-11.2(1) 


-12.9(1) 


-18.8(1) 


-23.8(1) 


<^e+i/e,PA/% 


-5.3(1) 


-8.6(1) 


-8.7(1) 


-9.2(1) 


-11.7(1) 


-14.1(1) 




-2.9(1) 


-5.3(1) 


-7.3(1) 


-9.0(1) 


-14.2(1) 


-18.5(1) 


'^^t+i/^,PA/% 


-2.8(1) 


-4.2(1) 


-4.8(1) 


-5.2(1) 


-7.1(1) 


-8.8(1) 


Srec/% 


-1.8(1) 


-2.7(1) 


-4.8(1) 


-6.3(1) 


-10.4(1) 


-13.6(1) 




-1.7(1) 


-1.6(1) 


-2.3(1) 


-2.5(1) 


-3.3(1) 


-3.9(1) 


A(2±iMw)/% 


±0.1 


±0.1 


±0.1 


±0.1 


±0.2 


±0.2 


A(Mt,.,z)/% 


-0.0 


0.1 


0.1 


0.2 


0.5 


0.7 



Table 3: Integrated lowest-order pp cross sections ctq for different ranges in pt,« and corre- 
sponding relative corrections 6, exact and in PA. The QED scale uncertainty is illustrated 
by A(M), as described in thet text. 



the plots as long as resonant W production is possible; only for transverse momenta pt,z 
much larger than Mw/2 differences become visible. 

The corrections to the high-transverse momentum tail are illustrated in Table ^ As 
expected, the PA again becomes worse with increasing p^^, but owing to the lower lu- 
minosity and the smaller cross section at the Tevatron in comparison to the LHC it will 
be extremely hard to see any effect of the enhanced electroweak corrections at high pT,i- 
Concerning the QED factorization scale dependence, which is again quantified by A(M), 
the same remarks apply that have already been made for the LHC case. The QED scale 
uncertainty turns out to be negligible. 

5 Conclusions 

We have calculated the complete set of electroweak 0{a) corrections to the Drell-Yan-like 
W-boson production at hadron colliders. Particular attention has been paid to issues of 
gauge invariance and the instability of the W bosons. All relevant formulae are listed in 
a form that facilitates their implementation in computer codes. Besides results for the 
full correction for off-shell W bosons, we have also presented an expansion of the virtual 
correction about the W-resonance pole, which is considerably simpler than the full result. 

Numerical results have been discussed at the parton level and for hadronic collisions 
at the LHC and at the Tevatron. The electroweak corrections significantly influence the 
transverse momentum and invariant-mass distributions of the decay leptons that are used 
in the determination of the W-boson mass. The pole approximation yields a good de- 
scription of the corrections to these observables. This result justifies, in particular, the 
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present practice at the Tevatron, where such an approximation is used in the W-mass 
measurement. However, in the domains of non-resonant W-boson production in these 
distributions, which are relevant for the measurement of the W-boson width or for the 
search of new-physics effects at the LHC, this approximation fails, rendering the complete 
correction important; the 0{a) corrections reduce the signal by several 10% for transverse 
lepton momenta with px,/ ^ 100 GeV. To further improve the analysis, in particular in 
the large-pT domain, more theoretical studies of electroweak higher-order effects as well 
as realistic experimental simulations are desirable. 

Note added 

While this paper was completed, another calculation of the electroweak 0{a) corrections 



to W production, including non-resonant contributions, has been presented in Ref. P2 
The numerical results mainly focus on observables measured at RHIC, but also include the 
transverse lepton momentum distribution at the Tevatron. However, no mass factorization 
has been performed in Ref. so that the corrections depend very sensitively on the light 
quark masses. A direct comparison with the numerical results presented in our paper is 
thus not possible. 

Acknowledgement 

We would like to thank U. Baur, M. Dittmar, W. HoUik, H. Spiesberger, and D. Wackeroth 
for helpful discussions. This work has been supported in part by the European Union under 
contract HPRN-CT-2000-00149. 

Appendix 

A Vertex and box corrections 
A.l Form factor for the Wff vertex 

The form factor used in (|2.12| ) for the vertex corrections is explicitly given by 



IQn I c4s4 si\ s 

+ - 24 + 2Q}si + 2Q),si) (2 + ^) 5o(0, 0, Mz) 

-4(1 + ^)5o(3, 0, Mw) - ^ f 1 + <LiM ) B,{s^ Mw, Mz) 



+ 



1 



2 



M' 

'i—iQf - Qf - 2Q/4)(Q/ - Qf + 2Q/'4) 



s 



+ 3(i-24-4g^g^,4) 



So(S,0,0) 
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+ 8Q}Bo{m}, 0,mf) + 8Q},Bo{m},, 0, m^O 

— 8QfQf'sCo{m'j, mj,, s,mf, m-y,mfi) 

+ 8Qf{Qf- Qf>)M^Co{ml 0, s, 0, rrif, Mw) 

-8Qf,iQf-Qf>)M^Coimj,,0,s,0,mf,,M^) 



+ 



8c: 



^'M^ + Ml + ^^^^ 



2 



Co(0,0,s,Mw,0,Mz) 



(Qf - Qf - 2Qfsl)iQf - Qf + 2Qf.s 



x[l + ^j sCo(0,0,s,0,Mz,0)j, 



(A.l) 



where the weak isospins I^ j and /w,/' the fermions are imphcitly taken to be Qf—Qf = 
21^ J = —21^ J,. Here and in the following the scalar integrals Bq, Cq, and Dq depend on 
their arguments as follows, 



Bo{pl,mo,mi) 



(27r/i) 



4:-D 



[g2 - mo + ie] [{q + pi^ - mj + ie] 
Co{pl,{p2-pif,pl,mo,mi,m2) = T^ f <fq 



[g2 -ml + ie] [(g + piY - mj + ie] [{q + P2Y -1712 + ie] ' 
Do{pl,{P2~Pif ,{P3-P2y, pI,pI,{P3- Pi f, mo, mi, 7712,7713) = ^ f d^q 



X 



[g2 - mp + ie] [(g + pi)^ - mf + ie] [{q + ^2)^ - ma + ie] [{q + pa)^ - + ie] 



(A.2) 



Explicit representations for the regular integrals can, e.g., be found in Ref. |]15|,|2^. The 
scalar integrals that involve mass-singular logarithms such as ln(m/) or on-shell singular- 
ities such as ln(s — + ie) are given by 



2 

So(m?,0,mi) = ^ + 2 + \n(^^, 



Bois,0,M. 



w 



A + 2 + In 



+ 



fM? 



w 



1 In 1 - 



Co{mi,m2, r, mi, m^, ma) 



In 

r L V mima 



m: 



7 



Co{ml,0,s,0,mi,Myj) 



4 ( m| 
^^nf^llnd 



\m1 



+ Li2 1 



^1/ 



— le 



2 



— le 



TT 
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Co{ml, 0, M^, m^, mi, Mw) 



1 i„f^£^)i„f^), 

V mi / VmiMw/ 



(A.3) 



where m^ <^ m^2 ^ M^. The quantity A = 2/(4 — — 7E + ln(47r) is the standard 
UV divergence for D — > 4 space-time dimensions, and fi is the (arbitrary) reference mass 
scale of dimensional regularization. 



A. 2 Box correction 

In order to reduce the Dirac structure of the box diagrams to the structure 



which appears in the lowest-order matrix element TWo) we used the identities 

[VdltnUJ-Uu] [Uut^u^^^Vi] = -tCo/2, 
[Vdl'^YY^-Uu] [u^^'J^'Julp^^^Vi] = 16Co, 



-4tco, 



which are valid in four space-time dimensions. 

The box correction factor 5box('S,t) introduced in ( |2.10| ) is explicitly given by 



^box(5,t) = £(s-m2,)| 



{9d9i +9u,9u) 



-{Boii, 0, 0) - Bois, Mw, Mz)) + 4Co(0, 0, s, Mw, 0, Mz 



u 



tCo(0, 0, t, 0, Mw, 0) + tCo(0, 0, t, 0, Mz, 0) 



+ 2sCo(0,0,s,Mw,0,Mz 
i{Ml + M| + s + 2if + 2(M2, + t)(M| + i)u 



(A.4) 



(A.5) 



Do(0,0,0,0,t, s,0,Mw,0,Mz) 



- 2{g^g-^ + gig-) [ 2Co(0, 0, S, Mw, 0, Mz) - nDo(0, 0, 0, 0, n, S, 0, Mw, 0, Mz) 
+ QdQi ^(5o(t,0,0)-5o(s,0,Mw)) +2Co(m2,0,s,0,m<,,Mw) 



+ 2Co(m2,0,s,0,m/,Mw) 



M^ + s + 2t 



(^sCo("^d,0,s,0,mrf, Mw) 



+ sCoirrii , 0, s, 0, m^, Mw) + tCo(0, 0, t, 0, Mw, 0) + tCo{m^, rrii , t, m^, m^, m^" 

2 , i\2^ 



t f 1 H j Do{ml, mf, 0, 0, t, s, m^, m^, m/, Mw 



2QuQi 



Co(mf , 0, s, 0, nil, Mw) + Co(m^, 0, s, 0, m„, M 



WJ 



uDo{m^, nil > 0) 0) "^7) "^/^ ^w) 



(A.6) 
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where u = — s — t, and are the left-handed couphngs of fermion / to the Z boson, 



- W,/ ■'w ^ r \ ^\ 



The singular Dq function is given by 

2 2 1 f 2 / "^1 \ 12/ "^2 

Z^oK,m,,0,0,r,.,m„m„m„Mw) = ^:p^j^|ln j + In 

with r = t, M < 0. 
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